Joule heating and thermoelectric properties in short single-walled carbon nanotubes: 

electron-phonon interaction effect 
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The electron-phonon interaction (EPI) effect in single-walled carbon nanotube is investigated by 
the nonequilibrium Green's lunction approach within the Born approximation. Special attention 
is paid to the EPI induced Joule heating phenomenon and the thermoelectric properties in both 
metallic armchair (10, 10) tube and semiconductor zigzag (10, 0) tube. For Joule heat in the metallic 
(10, 10) tube, the theoretical results for the breakdown bias voltage is quite comparable with the 
experimental value. It is found that the Joule heat can be greatly enhanced by increasing the 
chemical potential, while the role of the temperature is not so important for Joule heat. In zigzag 
(10, 0) tube, the Joule heat is smaller than the armchair tube, resulting from nonzero band gap in the 
electron band structure. For the electronic conductance G e and electron thermal conductance a e i, 
the EPI has important effect at higher temperature or higher chemical potential. Compared with 
ballistic transport, there is an opposite tendency for G e to decrease with increasing temperature after 
EPI is considered. This is due to the dominant effect of the electron phonon scattering mechanism in 
the electron transport in this situation. There is an interesting 'electron-drag' phenomenon for the 
phonon thermal conductance in case of low temperature and high chemical potential, where phonons 
are dragged by electrons from low temperature region into high temperature region through EPI 
effect. 

PACS numbers: 73.63.Fg, 72.10.Di, 72.20.Pa, 66.70.-f 
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I. INTRODUCTION 

The single-walled carbon nanotubes (SWCNT) were 
discovered in 1993^2, and have attracted a lot of re- 
search interest from both theoretical and experimental 
communities (for review, see e.g. Ref. 3-5). The SWCNT 
can be considered as a cylinder from rolling up of a two- 
dimensional graphene sheet in a particular direction, and 
it can be either metallic or semiconductor depending on 
the rolling direction. 6 Considering its good electronic per- 
formance, one promising application of the SWCNT is 
to design functional electronic nanodevice»2r— In these 
nanoscale electronic devices, the Joule heating due to 
the electron-phonon interaction (EPI) can be a disas- 
ter. This Joule heating is not only a waste of energy, 
but also can break down devices. There have been some 
studies on the in-air breakdown of metallic SWCNT by 
Joule heatingJ^r— We will visit this Joule heating topic 
in the SWCNT by the nonequilibrium Green's function 
(NEGF) approach in this paper. 

Another important physical property is the thermo- 
electricity, which is characterized by the figure of merit 
(ZT). The pursuit of high value of ZT has been in pro- 
cess for some decades, yet the ZT of most commonly 
used thermoelectric materials is still on the order of 
I 19.-22 j n recen t years, owning to the development of 
the nanotcchnology, some nano materials are found to 
exhibit high ZT, eg. the silicon nanowire with rough 
surface. 23 As the thermoelectricity is an important prop- 
erty, there have been large amount of studies focusing on 
this topici^r— In a recent experiment, it was proposed 
that the EPI should play an important role in the ther- 



moelectric properties ! 32 ' 33 However, to the best of our 
knowledge, there is still quite few theoretical works dis- 
cussing how the EPI can affect the thermoelectric prop- 
erties in the SWCNT. This is another topic of this paper. 

In this paper, we study the EPI effect in SWCNT 
by using the NEGF approach under the Born approx- 
imation. The metallic armchair SWCNT(10, 10) and 
semiconductor zigzag SWCNT(10, 0) are comparatively 
studied. We are interest in two phenomena: the Joule 
heating and thermoelectric properties affected by EPI. 
We find that in all tubes, the Joule heating is sensi- 
tive to the chemical potential applied, and only shows 
slightly dependence on the temperature. From the Joule 
heating we can obtain the in-air breakdown bias voltage 
Vbd- The theoretical result for Vbd agrees with the ex- 
perimental value. Compared with the metallic armchair 
SWCNT(10, 10), the semiconductor zigzag SWCNT(10, 
0) has smaller Joule heat, due to its nonzero band gap in 
the electron band structure. 

We also investigate the thermoelectric properties af- 
fected by the EPI. We find that the EPI effect is im- 
portant for the electronic conductance G e and electron 
thermal conductance a e i at higher temperature or higher 
chemical potential. The G e increases with increasing 
temperature in the ballistic transport. However, after the 
EPI is considered, we show that the G e exhibits decreas- 
ing behavior with the increase of temperature, resulting 
from the strong electron phonon scattering. For phonon 
thermal conductance, we find an interesting 'electron- 
drag' phenomenon at low temperature and high chemi- 
cal potential, where the phonon can be dragged by elec- 
tron from low temperature region into high temperature 



through the EPI effect. 

The present paper is organized as follows. In Sec. II, 
after a discussion of the electron and phonon Hamilto- 
nian, we present how to consider the EPI effect in NEGF 
approach within Born approximation. Subsec. Ill A is 
devoted to the results of Joule heat. The calculation 
method for Joule heat is discussed. Subsec. Ill B is for 
the results of thermoelectric properties. The numerical 
methods for different physical quantities are presented. 
The paper ends with a brief summary in Sec IV. 



II. THEORY AND FORMULA 



Phonon and electron Hamiltonian 
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FIG. 1: (Color online) Feynman diagram for EPI in the Born 
approximation, (a) is the phonon self-energy; (b) and (c) are 
the electron self-energy. 



is straightforward to show the following two symmetries 
of the M matrix: 



Km = M ml, 
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The lattice dynamic properties of the SWCNT are de- 
scribed by the empirical Brenner inter-atomic potential 34 
implemented in GULP—. GULP is applied to optimized 
the structure and generate the force constant matrix 
which is needed in the phonon Green's function. The 
electron Hamiltonian is obtained in the single 7r orbital 
tight-binding approximation scheme^: 
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where J = 3.0 eV and Ji = -6.0 eV/A are the onsite 
potential and hoping parameter. A and B are atoms in 
the two sublattices of SWCNT. ej g is the direction from 
atom (j, A) to its neighboring atom (j+S, B). The hoping 
parameter Jj is a negative value, which indicates that it 
is more difficult for electrons to hop between two atoms if 
the distance between them increases. After some algebra, 
we can obtain a general form for the EPI Hamiltonian of 
a system with N atoms: 
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where the summation of m is taken over the first-nearest- 
neighbors of atom I. And the third index, n, can either be 
atom I or m. In this general formula, the choice of atom 
I is arbitrary, not restricted to sublattice A anymore. It 



We note that Mp m is a vector with the upper script n as 
an atom index. It can be written in another form MP^*, 
where a = x,y,z. Or we will simply use notation Mp m 
with n as the n-th degree of freedom. 



B. EPI effect under the Born approximation 

In the Born approximation, the phonon/electron self- 
energy are expanded in terms of the EPI Hamilto- 
nian, and the expansion is truncated up to the sec- 
ond order^Zr— Fig. [T] are the three corresponding 
Feynman diagrams for the phonon/electron self-energy. 
The dot/solid lines corresponds to the phonon/electron 
Green's functions without EPI. If these non-perturbed 
Green's functions are replaced by full Green's functions 
with EPI, then the full Green's functions can be ob- 
tained self-consistently. This is the self-consistent Born 
approximation! 3 ^, It is a big challenge to apply the self- 
consistent Born approximation in a real system like 
SWCNT in this paper, since the self-consistent solution 
is very difficult to achieve for system with lots of degrees 
of freedom. There are also some other approximations to 
study the EPI effect. For example, in the lowest order ex- 
pansion approximation^Sr— it is assumed that the single- 
particle Green's functions and lead self-energies are inde- 
pendent of the energy. This assumption can simplify the 
integrals over energy, while retains the Pauli exclusion 
principle for fermionic particles. 

We apply the Born approximation so as to get valid 
results within reasonable calculation cost for SWCNT. 
Fig. [U (a) is the EPI induced phonon self-energy. The 
lesser and retarded versions of this self-energy will be of 
use: 
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where lit, is the self-energy due to the coupling of the 
system to the left lead. The right lead is analogous. 
In parallel, the electronic current from left lead is: 



Diagram (b) is the EPI induced electron self-energy. Its 
lesser and retarded versions are: 
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Diagram (c) is another EPI induced self-energy for the 
electron. Resulting from the loop in the diagram, the 
lesser version self-energy is zero. The retarded version 
for this self-energy is: 
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We can see this self-energy is actually a constant matrix. 
Its effect is simply to shift the electron energy level by a 
constant. 

After the EPI induced self-energy are obtained as 
shown above, the full Green's function of the system can 
be calculated from Dyson equations: 
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where II* 



n<, and £< = £<+ £<+£<, 
are the total phonon/electron self-energy. 

The thermal or electronic current from the leads 
into center can be formulated through the full Green's 
functional - — The phonon thermal current flowing from 
left lead into the system is: 
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The electron thermal current from left lead into the sys- 
tem is: 
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where /i£ is the chemical potential of the left lead. 



III. RESULTS AND DISCUSSION 

A. Joule heating 

In systems without EPI, the phonon thermal conduc- 
tance is purely due to the temperature difference of the 
left and right leads. In this situation, it can be shown that 
the phonon thermal current from two leads has the same 



ph 



magnitude while opposite in direction™; i.e I\ 
This simply means that the phonon thermal current runs 
out of one lead and jump into the other lead. We refer to 
this kind of thermal current as normal thermal current. 
In a system with EPI, the phonon thermal current can 
also be generated by applying temperature difference in 
left and right leads. However, besides this normal effect 
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FIG. 2: (Color online) SWCNT(10, 10). Joule heat calculated 
with different parameters of a. AT = 2aT. 

there is another method to induce the phonon thermal 
current. It can be driven by introducing bias voltage 
without temperature gradient across the system. In this 
situation, phonons are excited through scattering with 
moving electrons. This is the microscopic mechanism of 
the Joule heating. The Joule heat will transfer from the 
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system into two leads. The total Joule heat current is 
the collection of these heat currents: Jq = +1^ ■ We 
point out that this formula still holds even if there is nor- 
mal heat current due to temperature difference between 
two leads. Because the normal thermal current has the 



relation L 



ph 



R , as a result the normal heat current 



from left and right leads cancel with each other in the 
formula of Jq. This issue will be further discussed in the 
following. 



1. armchair SWCNT (10, 10) 

Fig. [5] shows the Joule heat Jq at different bias voltage 
V with T = 300 K and n = 0.4 eV. The Fermi energy 
e F = 0. The length of the SWCNT here is 1.5 nm with 
240 carbon atoms. In the experiment, it is advantageous 
to prepare SWCNT samples with length in micrometer 
scale; while the quantum mechanical NEGF theoretical 
approach prefers shorter SWCNT due to large memory 
required. So we focus on short SWCNT with length in 
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FIG. 3: (Color online) SWCNT(10,10). Joule heating V.S. V 
at different chemical potential (a), and at different tempera- 
ture (b). 

nanometer scale. The parameter a is the temperature 
difference ratio, i.e Tr — Tl = 2aT. It shows that the 
value of Joule heat does not depend on whether there is 
temperature difference between two leads or not. As we 
have mentioned in the above, this is because the normal 
phonon thermal transport due to temperature difference 
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FIG. 4: (Color online) SWCNT(10,10). Joule heat with dif- 
ferent chemical potential or temperature, (a). Joule heating 
V.S. T. (b). Joule heating V.S. p. 



of the left and right leads cancels with each other in the 
formula Jq — Jl + Jr. In the rest of the text, if we 
calculate the Joule heat, without mentioning it, we set 
a constant temperature for the system, i.e a — 0. In 
this case, the phonon thermal current calculated from 
Eq. ([20]) are purely the Joule heat. We note that in some 
figures of this manuscript, the quantity units are omit- 
ted. They are in terms of SI units: [T]=K, [/z]=eV, and 
[V]=V. ' 

Fig. [3] is the Joule heat V.S. bias voltage V. Panel 
(a) shows the result for a constant temperature 300K. 
If the chemical potential /i is zero or very small, the Jq 
increases very slowly with increasing bias voltage V, and 
the value of Jq is obvious only after V > 2 V. If the 
chemical potential /i is considerable large, eg. 1.4 eV, or 
2.1 eV, the Fermi surface goes into the second/third en- 
ergy level of the electron conduction band. As a result, 
the conductive electron density is greatly enhanced and 
much more phonons will be excited through EPI, lead- 
ing to rapid increase of Joule heat with increasing V. A 
large value of Joule heat can be achieved at bias volt- 
age V = 0.5 V. These theoretical results indicate that in 
short SWCNT with length L in the nanoscale, the elec- 
tron transport is almost in ballistic region even in the 
existence of EPI. So the electron carrier density need to 
be greatly enhanced so as to generate remarkably large 
Joule heat. Panel (b) shows the Joule heat V.S. V at 
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different temperatures. Only limited influence from the 
temperature can be found, since the energy level of the 
electron is in the order of eV. The temperature need to be 
in the order of 10 4 K to introduce considerable influence 
on the electron doping. 

Although the temperature has very small effect on the 
Joule heat, it is still interesting to investigate this effect 
in more detail. It carries some valuable information of 
temperature dependence for the EPI effect. Fig. @] (a) 
shows the Joule heat at different temperatures. If /i = 
and V = 1.0 V, the Joule heat is monotonically decreas- 
ing function of temperature. At T=1800 K, the value of 
Joule heat is reduced by 50% compared with that at 150 
K. At /i = 1.4 eV and V — 0.3 V, we can observe more 
abundant behaviors of the Joule heat in different tem- 
perature regions. Below 600 K, the Joule heat increases 
with increasing temperature; then reaches a maximum 
value around 600 K; and finally decreases with further 
increasing temperature. 

Fig. [5] compares the Joule heat per length in SWCNT 
(10, 10) with different length L at n = and T = 300 
K. As we can see the Joule heat per length is smaller in 
longer tubes, since the electron is almost in the ballistic 
transport region. However, if we check the Joule heat 
itself, we find that the Joule heat is larger in the longer 
tube than that in the shorter one. They are not exactly 
equal to each other. Especially in higher bias voltage 
region. The underlying mechanism is that electrons in 
longer tube spend more time in the system during its 
travel from one lead into the other. As a result, more 
phonon will be excited, leading to larger value of Joule 
heat. 
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the I-V curve. The experimental value of Vbd is about 
2.6 V for SWCNT samples with length 10 ± 5 nm, and 
diameter distributed in the range of [1.5, 2.5] nm. In 
our theoretical simulation, the SWCNT (10, 10) has the 
length as 3.0 nm and diameter 1.4 nm. They are not ex- 
actly in the same size of experimental samples. However 
we would like to do some qualitative comparison between 
our theoretical results and the experimental data. First 
of all, the in-air breakdown temperature for SWCNT is 
about Tbd = 873 K, which is due to Joule heating by 
applying bias voltage. This value of Tbd is obtained by 
the thermogravimetric analysis experiments! 44 ' 45 . In ad- 
dition, the breakdown of SWCNT is taken place in its 
middle region by Joule heating, and the temperature in 
this middle region can be obtained from the heat con- 
duction equation 1 ^: 

T = T Q + P'/g, (23) 

where Tq is the enviroment temperature, and g rs 0.17 
WK _1 m _1 is the net heat loss rate to the substrate per 
unit length. P' is the Joule heat per unit length. So the 
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FIG. 5: (Color online) SWCNT(10, 10). Joule haet per length 
in SWCNT with different lengths. 

Now we do some comparison of our theoretical results 
with experiments. In experiments, the metallic SWCNT 
can be put on insulating substrate, and applied a bias 
voltage on it After the bias voltage is increased to 
a certain value, the SWCNT will breakdown and this 
breakdown bias voltage (Vbd) can be recorded through 



FIG. 6: (Color online) SWCNT(10, 0). Joule heat V.S. 
bias voltage V with (a) different chemical potential, and 
(b)different temperature. 

breakdown Joule heat per unit length can be obtained 
from this formula: P' BD = (T BD - T ) * g m 80 Wm" 1 . 
Corresponding to this value in Fig. [SJ we can obtain the 
theoretical breakdown bias voltage Vbd ~ 2.5 V, which 
is quite comparable with the experiment. 
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2. zigzag SWCNT (10, 0) 

Fig. |5] is the Joule heat V.S. V in the zigzag 
SWCNT(10, 0). It has a similar behavior as the arm- 
chair SWCNT(10, 10) shown in Fig. O In case of small 
(j,, the Joule heat keeps small upto V = 1.0 V. If chem- 
ical potential is larger, the Joule heat increases quickly 
with the increase of bias voltage V. The temperature 
only plays a weak role on the Joule heat. Besides these 
similar behaviors, we can see some interesting difference 
between these two types of SWCNT. The SWCNT(10, 
0) has a band gap of about 1.2 eV in its electron band 
structure, while SWCNT(10, 10) is gapless. So the for- 
mer is semiconducting while the later is a metallic tube. 
As a result, in small n region the Joule heat in the zigzag 
SWCNT(10, 0) should be much smaller than that in the 
armchair SWCNT(10, 10). This is confirmed in our cal- 
culation shown in Fig. 0(a) and Fig.@](a) with /.i = and 
V = 1.0 V. In whole temperature range, the Joule heat 
in zigzag SWCNT(10, 0) is smaller than the armchair 
SWCNT(10, 10) at least by a factor of four in magni- 
tude. 

Fig. [7] (a) also shows that the Joule heat is almost 
zero at low temperatures in zigzag SWCNT, result- 
ing from the nonzero electron band gap. The Joule 
heat increases monotonically with increasing tempera- 
ture, which is quite different from the armchair SWCNT 
shown in Fig. [7J (a) with fi = and V = 1.0 V. As shown 
in panels (b) and (c), if the fi is a small nonzero value, 
the Joule heat shows a small decreasing behavior in low 
temperature region; then reaches a minimum value and 
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FIG. 7: (Color online) SWCNT(10, 0). Joule heat V.S. 
perature at different fixed chemical potential. 



increases with further increase of temperature. Panel (d) 
displays the result for large chemical potential, where the 
Joule heat decreases monotonically with the increase of 
temperature. 

Fig. [5] shows that the chemical potential also has im- 
portant effect on the Joule heat in the zigzag SWCNT. 
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FIG. 8: (Color online) SWCNT(10, 0). Joule heat V.S. chem- 
ical potential with different fixed chemical potential and bias 
voltage. 



If the chemical potential is high, a large Joule heat can 
be obtained with quite small bias voltage. 



B. thermoelectric properties 

The thermoelectric physical quantities can be calcu- 
lated by definition from the thermal and electronic cur- 
rents. This numerical calculation procedure is different 
from the ballistic transport where Landauer formula can 
be used for the thermal and electronic current. Actually, 
the expression in Eq. (|2T)1) - (|2"2"|) can be reformulated into 
an effective Landauer formula after an effective trans- 
mission function is introduced. 39 However, this effective 
transmission function is quite different from its counter- 
part in the ballistic transport, because it depends on 
temperature (for phonon and electron) and also chemical 
potential (for electron). As a result, we need to do calcu- 
lation numerically for the first derivative of the thermal 
or electronic currents with respect to the temperature or 
chemical potential. These first derivatives of the ther- 
mal/electronic currents are required in the calculation of 
different thermoelectric physical quantities. In ballistic 
transport, they can be simply obtained by taking first 
derivatives of the Bose-Einstein or Fermi-Dirac distribu- 
tion functions. 

The Taylor expansion of the thermal and electronic 
currents in terms of AT = Tr — Tl and Afj, = fJ,R — fJ>L 
are 
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FIG. 9: (Color online) SWCNT(10, 10). Thermoelectric prop- 
erties calculated with different parameters (a, 0), a and /3 are 
the ratio of difference for the temperature and chemical po- 
tential, respectively. 



where the partial differential coefficient can be obtained 
numerically, for example, 



dI ph p h (T, AT, fi, Afi = 0) 



dT 



AT 



(27) 



Using these numerically obtained coefficients, the ther- 
moelectric quantities can be calculated by definition. The 
phonon thermal conductance is: 
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(28) 



The electronic conductance is: 
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FIG. 10: (Color online) SWCNT(I0, 10). The EPI effect on 
the thermoelectric properties at fixed temperature. 
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1. armchair SWCNT (10, 10) 

As mentioned above, the partial differential coefficients 
need to be calculated numerically. A small temperature 
difference AT, and chemical potential difference A/i be- 
tween left and right leads are introduced. The temper- 
atures in the left and right leads are Tl — (1 — a)T 
and Tr = (1 + a)T ', where T is the average tempera- 
ture. The chemical potentials are (j-l = (1 — P)n and 
lift i = (l + /3)/i with n as the averaged chemical potential. 
So AT = 2oT, and Afi = 2/3/x. Smaller value of (a, 0) 
has the benefit of approaching the analytical partial dif- 
ferential coefficients numerically, yet too small value of 
(a, j3) may cause larger numerical error. So a proper cho- 
sen value for (a, /?) is important. Fig. [SJ compares two 
results for the thermoelectric quantities, with respect to 
(a,P) = (-0.1,-0.1) and (-0.05, -0.05). The chemical 
potential /i = 0.4 eV. It turns out that these two curves 
are almost indistinguishable. So (a, (3) = (—0.1,-0.1) is 
a good choice and will be used in following calculation 
without special mention. 

Fig. [TU] illustrates the EPI effect on different thermo- 
electric properties as a function of temperature, with 
chemical potential fx = 2.1 eV. Compared with the bal- 
listic transport, these physical quantities are generally 
depressed by the EPI. In low temperature region, the 
EPI only has small effect due to few energy and charge 
carriers, especially quite few phonons. With increasing 
temperature, more electrons will be scattered by phonons 
through the EPI, resulting in further reduction of these 
physical quantities. Panel (a) shows that the electronic 
conductance increases with the increase of temperature 
in ballistic transport, but it decreases after the EPI is 
considered. This interesting effect will be further inves- 
tigated in more detail in the following. Panel (c) and (d) 
exhibit that the Seebeck coefficient and ZT are slightly 
affected by the EPI around room temperature. So in the 
following text we are not going to spend much discus- 
sion on these two quantities, instead we will focus on the 
electronic conductance G e , electron thermal conductance 
a e i, and the phonon thermal conductance cr p h- 

In Fig. [TTJ we do further studies on how EPI affect 
the G e with different chemical potential. There are two 
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FIG. 11: (Color online) SWCNT(10, 10). The electronic con- 
ductance V.S. temperature, with fixed chemical potential. 

competing mechanisms on the G e with increasing tem- 
perature. On the one hand, more electrons are involved 
in transporting of electronic charges with the increase 
of temperature. It makes positive contribution to G e . 
On the other hand, more phonons will also be excited 
at higher temperature, leading to stronger scattering of 
electron by phonon. This mechanism has negative con- 
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FIG. 12: (Color online) SWCNT(10, 10). The electronic con- 
ductance V.S. chemical potential, at constant temperatures. 

tribution to G e . These two competing mechanisms result 
in rich dynamics of the G e . Panel (a) is for small chemical 
potential /i = 0.4 eV, which corresponds to low electron 
doping. In this situation, the electron number is too small 
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FIG. 13: (Color online) SWCNT(10, 10). The electron ther- 
mal conductance V.S. chemical potential, with different fixed 
temperatures. 



to have significant electron phonon scattering. As a re- 
sult, the EPI only has small effect and the positive mech- 
anism is more important than the negative mechanism. 
Consequently, G e is slightly smaller than that in the bal- 
listic transport, and keeps on increasing with increasing 
temperature. Panel (b) is for \i = 1.4 eV, where the EPI 
effect becomes more important and the negative mecha- 
nism begins to compete with the positive mechanism. As 
a result, G e exhibits obvious deviation from the ballistic 
result, and the increasing velocity is slowed down. Panel 
(c) is for large chemical potential /i = 2.1 eV. This is 
a high electron doping with Fermi level shifted into the 
third electron conduction band. In this case, the EPI is 
very important, and the negative mechanism dominants 
the electron transport property. So G e shows decreasing 
behavior with the increase of temperature. 

Fig. [12] illustrates the EPI effect on the behavior of G e 
versus /i. At all studied temperatures, the EPI has im- 
portant effect only in high chemical potential region. The 
EPI effect on the G e is almost ignorable in case of small 
chemical potential due to small number of electron. A 
very similar result is found for the electron thermal con- 
ductance a e i as shown in Fig. [T3l where the EPI effect is 
important only if the chemical potential shifts the Fermi 
level into higher electron conduction band. 

Fig. 03] shows the EPI effect on the phonon thermal 
conductance with different fixed chemical potential. As 
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FIG. 14: (Color online) SWCNT(10, 10). The phonon ther- 
mal conductance V.S. temperature, with different chemical 
potential. 



well known, a p h is independent of the chemical potential 
in the ballistic transport. However, by considering the 
EPI effect, the chemical potential can affect the phonon 
thermal transport through electron phonon scattering. 
This figure shows that a p h decreases obviously with in- 
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potential (J, = 2.1 eV, a p h is almost zero in low tem- 
perature region, which means phonons are almost com- 
pletely scattered by electrons. The reason is that only 
few phonons are excited at low temperature, yet there 
are huge number of conduction electrons at high chemical 
potential. As a result, all phonons are scattered by elec- 
trons, leading to small or zero phonon thermal conduc- 
tance. At this point it is interesting enough to imagine 
what will happen if the chemical potential is further in- 
creased. We can anticipate that the phonon thermal con- 
ductance may become negative, which means the phonon 
thermal energy transfers from low temperature region 
into high temperature region with the help of the electron 
phonon scattering. Indeed, this anticipation is confirmed 
by our calculation shown in Fig.QjjJ Panel (a) shows that 
the a p h at 150 K decreases with increasing chemical po- 
tential, and becomes negative after /i > 2 eV. It is worth 
noting that, this effect is a kind of 'electron-drag' effect 
and has nothing to do with the Joule heat, because in 
this calculation only very small bias voltage is applied 
V = A/i/e = 2/3 /i/e. The Joule heat under such small 
bias voltage is almost zero as can be seen from Fig. [3] 
This effect has something similar to the phonon-drag ef- 
fect, which describes charge carriers dragged by phonon 
from hot to cold regions via momentum transfer 4£ To 
see this 'electron-drag' effect on phonon thermal conduc- 
tance, high density of electron and low density of phonon 
are required; i.e high chemical potential and low temper- 
ature. Under these conditions, the SWCNT can serve 
as a heat pump to transfer phonon thermal energy from 
low temperature region into high temperature region by 
applying bias voltage. 
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FIG. 16: (Color online) SWCNT(10, 0). The thermoelectric 
properties calculated with different parameters (a, j3). 



FIG. 15: (Color online) SWCNT(10, 10). The phonon ther- 
mal conductance V.S. chemical potential, at constant temper- 
atures. 

creasing chemical potential. Especially for high chemical 



2. zigzag SWCNT (10, 0) 

The rest of the text is devoted to the thermoelectric 
properties in zigzag SWCNT(10, 0). Fig. [T51 supports 
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FIG. 17: (Color online) SWCNT(10, 0). The electronic con- 
ductance V.S. temperature, with different fixed chemical po- 
tential. 

that parameters (a,/3) = (—0.1,-0.1) are also suitable 
choice in zigzag SWCNT(10, 0). The chemical potential 
/i = 0.3 cV in the figure. We note that the Seebeck coeffi- 
cient in zigzag SWCNT is 2-3 orders of magnitude larger 
than the armchair SWCNT, resulting from the band gap 
in the electron band structured For the same reason, 



zigzag SWCNT. As a result, the figure of merit ZT is 
considerably larger in zigzag SWCNT. 

Fig. [17] is the electronic conductance G e in 
zigzag SWCNT(10,0) with different chemical poten- 
tial. Compare this figure with the result of armchair 
SWCNT(10,10) shown in Fig. QU some similar and dif- 
ferent behaviors can be found. In both figures, the G e in 
ballistic transport will increase with increasing temper- 
ature. If the EPI is taken into account, the results de- 
pend on the value of chemical potential. For low chemical 
potential/i = 0.3 eV, G e only deviates slightly from bal- 
listic result in low temperature region; and this deviation 
increases with increasing temperature, as more phonons 
are excited at higher temperature leading to stronger 
electron-phonon scattering. However, G e keeps on in- 
creasing as the temperature increases, which is same as 
ballistic transport. For high chemical potential /i = 0.9, 
1.5 eV, the effect of EPI becomes more and more impor- 
tant, because the possibility of electron-phonon scatter- 
ing increases as the chemical potential increases. So the 
G e is much smaller than the ballistic result, especially in 
the high temperature region, and G e will decrease with 
increasing temperature which is opposite from the bal- 
listic case. The difference between Fig. [17] and Fig. QT] 
is quite obvious. At low chemical potential, the G e in 
SWCNT(10,0) is almost zero at low temperature, and 
is much smaller than that in SWCNT(10,10) in whole 
temperature range. This is because of the band gap in 
the electron band structure in the zigzag SWCNT(10,0), 
so the electronic conductance is mainly contributed by 




FIG. 18: (Color online) SWCNT(10, 0). The electronic con- 
ductance V.S. chemical potential, at constant temperatures. 

the electron thermal conductance is much smaller in the 



the doping electrons in low temperature region. Due to 
the same reason, the EPI effect is more obvious in the 
SWCNT(10,0). For example, in panel (c), the reduction 
of G e can be as much as 50% by the EPI, which is quite 
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FIG. 20: (Color online) SWCNT(10, 0). The phonon thermal 
conductance V.S. temperature, with fixed chemical potential. 



ductance, since the thermal current and electron charge 
current are carried by the same carriers (electron). For 
all three considered chemical potential, the EPI has lim- 
ited effect in low temperature region and becomes very 
important in high temperature region. This result brings 
some information that, the EPI is more sensitive to the 
temperature in semiconductor zigzag tubes than that in 
the armchair tubes. 

The phonon thermal conductance versus temperature 
is shown in Fig. [20] It is quite similar as the arm- 
chair tubes, because the difference between these two 
types of tubes is in the nonzero band gap of the zigzag 
SWCNT(10, 0). However this band gap has no direct 
effect on the phonon thermal transport. The 'electron- 
drag' effect for the zigzag SWCNT(10, 0) can also be seen 
from Fig. [gl] (a) where T = 150 K. After fi > 1.6 eV, the 
phonons transport from low temperature region into high 
temperature region, driven by the dragging force of elec- 
trons through electron phonon scattering. 



larger than that in the armchair tube. 

Fig. [18] shows the importance of temperature and 
chemical potential in the SWCNT(10,0). In low chemical 
potential region, the EPI is not so important; while the it 
becomes very important at high chemical potential. At 
1000K, the G e is reduced by almost one half. 

Fig. Q1J] is for the electron thermal conductance versus 
temperature at different chemical potential. This figure 
is more or less similar as Fig. [T7] for the electronic con- 
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FIG. 21: (Color online) SWCNT(10, 0). The the phonon ther- 
mal conductance V.S. chemical potential, at constant temper- 
atures. 



IV. CONCLUSION 

To conclude, the EPI effect in SWCNT is investigated 
by using the NEGF approach within the Born approx- 
imation. The Joule heat and thermoelectric properties 
in both armchair SWCNT(10, 10) and semiconductor 
zigzag SWCNT(10, 0) are comparatively studied. It was 
found that the Joule heat is more sensitive to the chemi- 
cal potential (gate voltage) than the temperature. With 
higher chemical potential, the Joule heat can reach a very 
large value by applying small bias voltage. The break- 
down bias voltage for armchair SWCNT(10, 10) is esti- 
mated to be about 2.5 V, which is quite comparable with 
the experimental value 2.6 V. The Joule heat in zigzag 
SWCNT(10, 0) is considerably smaller than the armchair 
SWCNT(10, 10), resulting from the nonzero band gap in 
the electron band structure. 

For the thermoelectric physical quantities, it was found 
that the EPI plays an important role at higher tempera- 
ture or higher chemical potential for the G e and a e i . At 
high chemical potential, the G e decreases with increasing 
temperature, which is opposite of the ballistic transport. 
This is the result of very strong electron phonon scatter- 
ing. The phonon thermal conductance exhibits an inter- 
esting 'electron-drag' phenomenon in case of low temper- 
ature and high chemical potential. In this situation, the 
phonons can be dragged by electrons from low temper- 
ature region into high temperature region with the help 
of EPI. 
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